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ABSTRACT 

The rate evolution of subluminous Type la Supernovae is presented using data from the Supernova 
Legacy Survey. This sub-sample represents the faint and rapidly-declining light-curves of the observed 
supernova la (SN la) population here defined by low stretch values (s < 0.8). Up to redshift z = 0.6, 
we find 18 photometrically-identified subluminous SNe la, of which six have spectroscopic redshift 
(and three are spectroscopically-confirmed SNe la). The evolution of the subluminous volumetric 
rate is constant or slightly decreasing with redshift, in contrast to the increasing SN la rate found 
for the normal stretch population, although a rising behaviour is not conclusively ruled out. The 
subluminous sample is mainly found in early-type galaxies with little or no star formation, so that 
the rate evolution is consistent with a galactic mass dependent behavior: r(z) —Ax M g , with 



A = (1.1 ± 0.3) x 10" 14 SNeyr 
Subject headings: supernovae: general 



1. INTRODUCTION 

While it is commonly agreed that Type la Super- 
novae (SNe la) are thermonuclear disruptions of mass 
accreting C-0 white d warfs (WDs) in a binary system 
(|Hovle &i Fowler|[i960D . the physics of the explosion and 
the nature of the companion are still under discussion. 
The study of the properties and environments of SNe la 
can help us solve the progenitor question. Successful pro- 
genitor and explosion models need to explain the variety 
in light-curve shape and spectra of SNe la. Recent multi- 
dimensional simulation s of the explosion (iGamezo et al 



2005; Liv ne et al.ll2lM IKuhlen et aLl |t006; Kasc n et al 



20091) show the asymmetric character of the delayed det- 
onation and succeed to explain the scatter in the width- 
luminosity relation of the normal SNe la, although com- 
plications for the extreme SNe la still exist. In particu- 
lar, subluminous SNe la - a group of objects considerably 
fainter at peak (up to 2 mag), with faster light-curves, 
redder at early phases, with distinct spectral character- 
istics such as Ti II and enhanced Si II - pose a challenge 
to any successful progenitor and explosion theory. 

gonzalez@astro.utoronto.ca 



The prototyp ical example of a sublu minous SN la is 
SN1991bg (e.g. IFilippenko et al.lll99"l but many other 
examples have been discovered and studies of their prop- 
erties relative to the bulk sample have been under- 
taken (iGarnavich et all 120041: iTaubenberger et al.ll2Ci08t 
iKasliwal et all 120081: iHicken et al.ll2009f ). Subluminous 
SNe la are predominantly found in galaxies with older 
stell ar population s, such as ellipticals and early-type spi- 
rals (jHowcll 2001), in contrast to the preference for the 
late-type, star f orming hosts favoured by bright, slow de- 
clining SNe la (|Hamuv et all 119961 [2000 ). Subluminous 
SNe la also occur exclusively in massive galaxies whereas 
the normal samp le spans a wider host stellar mass range 
(jNeill et al.ll2009D . These observations set constraints on 
the delay-time - the time between the formation of the 
binary system and subsequent SN explosion - ranging 
from < 1 Gyr for SNe occurring in star-forming regions 
to several Gyrs for SNe in quiescent environments. They 
could as well hint at differences in the metallicity abun- 
dance of the progenitors, du e to the mass-metallicity re- 
lation (jTremonti et al.l [20041 ). i. e., subluminous SNe la 
happen in more metal-rich environments. 

Based on a sample from the Lick Observatory Super- 
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nova Search (LOSS) and the Beijing As tronomical Ob- 
servat ory Supernova Survey (BAOSS), iLi et all (|200lL 
[20101 found that 17.9±£J% of all local SNe la are 
1991bg-like, a possible overestimation as these surveys 
were host-targeted and could have systematically sam- 
pled brighter and more massive galaxies. The current 
observed subluminous sample remains a small fraction 
of the overall SN la population, although new and re- 
cent surveys are actively looking for them (jHicken et al.l 
2009). At high redshift, they are challenging to iden- 
tify spectroscopically, and the current high-z SN surveys 
usually preferentially target normal SNe la for cosmolog- 
ical purposes. As a result, no 91bg-like SNe have been 
located at z > 0.1. For example, both, the Supernova 
Legacy Survey (SNLS) and the ESSENCE survey, report 
no spectra of 91bg-like objects at z > 0.1 although this 
is consistent with thei r selection effects (Brondc r et al.l 
2008t iFolev etHI 120091 ) . 

This raises an obvious question: is this non detection 
simply due to the difficulty of detecting and classifying 
these objects, or is the relative frequency actually lower 
at high redshift because there has not been e nough time 
for them to explode as SNe la (|Howellll200"TI )? Problems 
in the detection and classification arise from their intrin- 
sic faintness, their rapidly evolving light-curves (which 
cause them to spend less time above the detection thresh- 
old for the same brightness) , and their tendency to occur 
in brighter galaxies where the low contrast between the 
SN and host makes spectroscopy difficult. In this paper, 
we look for the fastest (and therefore faintest) SNe la 
at z > 0.1 in th e Supernova Legacy Survey (SNLS) 
(|Astier et al.l 12006). Even without spectroscopic follow- 
up, we can use the excellent multi-band light-curves of 
the SNLS to look for subluminous SN la candidates by 
fitting subluminous SN la LC templates to the photo- 
metric data. 

Beyond a simple detection, even more enlightening 
would be a measurement of the evolution in the volumet- 
ric rate of these objects. As it is a convolution of the star- 
formation history and the delay-time distribution, the 
SN la rate, in particular for sub-samples of SNe la, can 
test different progenitor scenarios. Rat es for SNe la have 
been measured by several gr o ups (e.g. IPain et al.l I 1996L 



2001 ICappellaro et al.l [19991: IDahlen et all 12004 120081: 
Neill et al.l 120061: iDildav et al.l 120081 120101: iPerrett et al l 
2010aHLi et al.ll2010bl) . Combining the results one finds 
an increase with redshift up to z ~ 1.0, although with 
large spread among surveys. If the rate evolution for 
the individual SN la populations differ, their delay-times 
must consequently vary and imply different conditions for 
the progenitors. As we probe higher redshifts, the SN la 
environments will differ from the local sample, and will 
be reflected in the rates. If the fraction of high-mass (and 
high-metallicity) hosts was lower in the past, we would 
expect a similar behaviour for the rate of subluminous 
SNe la. 

By constraining the delay-time distribution, we seek to 
improve our understanding of the progenitors of sublu- 
minous SNe la. One can attempt to explain their progen- 
itors in the SN la progenitor framework, where two main 
scenarios have been envisaged: the single-degenerate 
(SD) model, in which a non-degenerate companion do- 
nates H/He-rich material to a WD near the Chan- 



drase khar mass (|Whelan fe Ibenl Il973t iNomoto et al.l 
19841) . and the double-degenerate scenario of tw o WDs 
coalescing (jlben fe Tutukovlll984lWebbinklll984l) . In the 
SD scenario the delay-time is set by the age of the donor 
while in the DD model it depends on the age of the sec- 
ondary and the merging time of the two WDs through 
gravitational wave radiation. The merging of two WDs 
has recently shown to be a viable mechanism, under cer- 
tain conditions, for a successful subluminous explosion 
(jPakmor et al.H2010fl . 

There is also a variety of independent mechanisms for 
subluminous explosions that explain the faintness by re- 
quiring a smaller burned mass of 56 Ni. Instead of a com- 
mon explosion mod el for the whole ran ge of SNe la, the 
delayed detonation (jMazzali et al.ll200"7j) and other mech- 
anisms like pure deflagrations or sub-Chandra sekhar 
CO or O-Ne-Mg WD off-center detonations (e.g. ILivnel 
19951: IWooslev fc Weaven fl99l iFilippenko et al.l 11992 



Isern et aTT 19911 ) provide less lum inosity although with 



some observational discrepancies (lHoeflich fc Khokhlovl 
1996). A failed neutron star model, in which C-0 is 
accreted rapidly and ignite d on the surface of th e WD 
also leads to a faint SN la (No moto fc Ibenl I1985D . Re- 
cently proposed models of WD collisions in dense en- 
viron ments could potentially lead to sublumino us explo- 
sions ([Raskin et al.ll2009at iRosswog et al.ll2009D . Finally, 
a theoretically-motivated group of unusu al objects, the 
so call ed SNe .la ("dot la") introduced bv lBildsten et al.l 
■ 2007;. are the final helium shell flashes of AM CVn bi- 
naries. These are too dim and evolving too rapidly to be 
classical sublumino us SNe la, but possible candidates are 
start ing to appear (Poznans ki et al.ll2010l : Kasliwa l et al.l 

With a detection of subluminous SNe la at z > 0.1, 
the measurement of the rate evolution and study of 
their host environments, we can start setting constraints 
on the evolution of these objects and their proper- 
ties, their delay- times and progenitor s. This paper is 
a companion to IPerrett et al.l (|2010a[) . who have car- 
ried out a rate study of normal SN la sub-samples us- 
ing the same SNLS data sample. We complement that 
study for the faint and short-lived SN la population 
up to z = 0.6. We use "stret ch" as a LC-shape pa- 
rameter dPerlmutter et"aLl ll997t iGoldhaber et all 120011: 
IConlev et alT 2008) and define hereafter subluminous ob- 
jects as SNe la with s < 0.8, inspired by the behav ior o f 
the color-stretch and magnitude-stretch relations f ^2.1[) . 
This value is higher than the ones inferred from othe r 
studies (e.g. iTaubenberger et al.l 120081 : ILi et al.l [^OlOa). 
as will be shown later. We emphasize therefore that our 
sample is less extreme and call it hereafter "low-stretch" 
(low-s) sample. 

A brief outline of the paper follows. In §21 we inves- 
tigate the nearby training sample of low-s SNe la and 
find the necessary LC relations for creating a template, 
^highlights the observations and candidate selection of 
the SNLS. With this, we measure the rate as a function 
of redshift in [21 taking into account systematic effects 
from possible errors in the estimated LC parameters and 
contamination from core-collapse supernovae (CC-SNe). 
The discussion of the rate measurements together with 
the investigation of the host galaxy properties of the can- 
didates follow in S}5] and we conclude in |j6j We use a flat 
cosmology with (fi TO , ti A ) = (0.27,0.73) and h = 0.7. 
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2. LOCAL RELATIONS FOR LOW-STRETCH SNE IA 



As shown by 
ITaubenberger et all (120081) 
and IHicken et all ([20090 , 



Garnavich et al.l ([2004 
IKasliwal et all ([2008 
subluminous SNe la 



have 

clear photometric characteristics: red colors, a faint 
peak brightness, and a steeper relationship between LC 
shape and peak magnitude. In this section we derive 
new parameter relations based on a training sample of 
subluminous objects at low redshift. These relations will 
be of great importance to define proper low-s SN la LC 
templates and to correctly calculate the SNLS survey 
efficiency. 

We use a local training sample from the literature (the 
low- z Constitution set and low-s, Table [T]) that spans 
the stretch range, 0.45 < s < 1.30, so that the behavior 
of low-s SNe can be compared with the normal popu- 
lation. Adequate light-curve phase coverage is ensured 
by requiring observations in at least two filters between 
-15 and 8 and between 5 and 25 effective days (corrected 
for stretch and redshift) after maximum brightness in 
B-band. We also include a sample at higher redshift 
from the SNLS consisting of spectroscopically-confirmcd 
SNe la over 0.3 < z < 0.6, which provides wider wave- 
length coverage compared to the low-z sample. 

We perform SN la light curve fits to the entire sam- 
ple with a modifie d version of the SiFTO algorithm 
([Conlev et al.ll2008h . SiFTO manipulates spectral energy 
distribution (SED) templates to model SN la LCs. It has 
a wavelength-stretch dependence derived from strictly 
normal SNe la, which was removed for the modified ver- 
sion used here. Although SiFTO was built with a train- 
ing sample that excluded peculiar objects like SN1991bg, 
we obtain very good fits to subluminous data down to 
s ~ 0.65 and reasonable fits to lower stretch SNe (see 
Figure [1} . Surprisingly, this template performs better 
than those built purely from subluminous objects like 
SN1991bg, which may have extreme features not suit- 
able for all s < 0.8 objects. The success of the fitting 
is partly because SiFTO fits each band independently 
without imposing any color-stretch relation, and so the 
unusual colors of subluminous SNe la do not pose a prob- 
lem. The algorithm provides the LC parameters such as 
stretch and peak magnitudes to get colors at maximum 
B-band brightness with their respective errors, as well as 
quality-of-fit parameters. 

We use a reduced overall goodness-of-fit xt < 5 cu t f° r 
normal and a higher xi < H f° r low-s objects to account 
for the poorer fits of the more extreme objects. The final 
training set consists of: 207 low-z SNe la, of which 44 
are low-s (see Table [IJ with s down to s = 0.45, and 106 
SNLS spectroscopic SNe la of which 3 are low-s. A list of 
all culls is shown in Table [3] We note th at not all of the 
low-s SNe la listed bv INeill et afl (f2009h (Table 1) pass 
our selection criteria. Also, their stretch values for low-s 
objects are in average lower than ours. This is due to the 
modified version of SiFTO we use, which assumes a com- 
mon single stretch independent of wavelength and which 
has been more extensively tested for low-s SNe la. Ad- 
ditionally, despite our light curve coverage criteria, one 
low-s candidate (SN2001da) does not have many data 
points and results in very diffe rent stretch (s = 0.89) 
when fit through the method of INeill et all (|2009f> . We 
take this supernova out of our analysis. 



In order to develop tools to photometrically identify 
low-s SNe la and distinguish them from other transients, 
we also include a nearby and a high- z SNLS sample of 
spectroscopically-identified non-type la supernovae. 14 
nearby and 21 SNLS CC-SNe are properly fit to a SN la 
SiFTO template, and 10 of those are assigned a s < 0.8 
stretch (6 of which are low- 2, see Table [5]). This set 
of contaminants will give us an idea of the systematic 
uncertainties (more in £|4.3[) . 

Our limiting definition of low-s objects at s = 0.8 is in- 
spired by the color-stretch and magnitude-color-stretch 
relations and is different from the "subluminous" def- 
inition studied in other surveys. We ca n compare it 
to the SALT2 LC-fitt er dGuv et all 12007ft using similar 
equations to Eq. 6 in lGuv et al.l ([2010D . obtaining that 
s = 0.8 for our modified SiFTO version corresponds to 
A~ i ~ —1.85. By directly m atching subluminous SNe 
of ITaubenberger et al.l ([2008ft and fitting a linear rela- 
tion between their Ami^(B) parameter (the amount of 
-B-band magnitudes the SN has declined 15 days after 
maximum) and our SiFTO stretch, we obtain a corre- 
sponding value of Ami^(B) ~ 1.54. This value is lower 
than the one they find to separate their sample (~ 1.70, 
which would correspond to s ~ 0.72) and confirms again 
that our sample is not as extreme. 
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Fig. 1.— Example SiFTO LC fits of low-s SN2005am 
IIHicken et alj 120091 ; [Li et al.ll2006t) with stretch 0.71 and overall 
X 2 V = 0.6 and SN1999by (lower) with stretch 0.61 and xl = 8.5. 
The fit is done for all points between -15 and 35 days in restframe. 



2.1. Color-stretch relations 
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Modeling the colors of SNe la is fundamental for the 
LC fit and candidate selection of SJ21 as well as for the 
study of the recovery efficiency needed for the rate cal- 
culation of [jU The challenge in modeling SN la colors 
resides in the difficulty of disentangling intrinsic and ex- 
trinsic colors, the latter coming m ainly from reddenin g 
of the host galaxy. As explained in lConlev et al. (2008), 
the SiFTO method does not impose a single color model 
across all filters, which would probably not encompass 
the whole range of SN la variability; rather the SED is 
adjusted to match the observed colors in each filter. 

We first remove highly red dened objects by the fore- 
ground Milky Way using the iSchlegel et all (| 19981) dust 
maps and a cut of E(B — V^)mw < 0.15. To get an idea 
of the intrinsic color-stretch relation of the objects, we 
take the transients that are hosted in less extincted pas- 
sive galaxies, i.e. in galaxies catalogued as E/SO for the 
nearby sample from the NED databasjB and in galax- 
ies wit h no star formation for the SNLS (Su llivan et al.l 
2006b). A color similar to B — V at B-band maximum 
brightness, c, as a function of stretch for these objects is 
shown in Figure [5] (filled circles and squares) . 

We create a "low-extinction sample" by taking all SNe, 
regardless of their hosts, with a color below the upper 
limit from the maximum found for the SNe in only qui- 
escent environments, c < 0.37. This is only done for 
normal-s SNe la, i. e. s > 0.8, which have a relatively 
flat color-stretch relation. Low-stretch SNe la become 
redder at a steeper rate with decreasing stretch and are 
hosted in less dusty environments, so we do not use such 
a maximum color cut, although we remove known ex- 
tremely extincted objects (e.g. SN1986G, SN2006br and 
SN2007ax). The fit to this larger sample (open and filled 
circles and squares) in Figure [5] is well represented by a 
piecewise linear model with a break at s — 0.77 (solid 
line): 



c(s) 



Co + 7i (s - 1) 



(s > 0.77) 



Co + 72 (a - 1) + C match (s < 0.77) 



(1) 



The fit parameters are Co = 0.08 ± 0.02, 7l = -0.06 ± 
0.02, 72 = —2.79 ± 0.07 and the transition stretch value 
s = 0.77 ± 0.01. Cmatch = —0.23(7i — 72) is a correction 
factor to unite both lines at s = 0.77. The transition 
value for stretch (s = 0.77), together with the one found 
for the magnitude-color-stretch relation in next section 
(s = 0.82), argues for our chosen limiting value of s ~ 0.8 
between low-s and normal-s. The piecewise linear model 
has the best fit (reduced xl 01 23.4) among other mod- 
els (e. g. exponential with xt = 25.0 or power law with 
xt = 25.9). The figure shows that low-stretch SNe la are 
not only redder but their color-stretch relation is much 
steeper than for normal SNe la. Errors in stretch and 
color were propagated in the fits. Allowing a fit to all 
objects, also extincted ones, provides consistent results 
(with a transition stretch at s = 0.83) and is also shown 
in Figure [2] It is worth to mention that if we force 
a piecewise linear model with a fixed lower transition 
value of s = 0.7, more similar to characteristic "sub- 

1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration 



luminous" surveys (jTaubenberger et ail 12008). the fit is 
slightly worse (xt = 24.6) and the slopes steeper. Due to 
the large scatter of colors at low-s, a range of consistent 



models (with X < X™ 
values of 0.69 - 0.87. 



2) is obtained for transition 




0.8 
stretch 

Fig. 2. — Color as a function of stretch for the local SN la sample 
(circles) and 0.3 < z < 0.6 SNLS sample (squares). The filled 
points represent SNe in passive galaxies and the fit to exclusively 
those objects is the dashed line, a two-piece linear model changing 
at s ~ 0.8. The "low-extinction" sample (open and filled black 
points) represent all objects with c < 0.37 regardless of the host 
galaxy. The solid line is the best two-piece linear fit to this sample 
and the dotted lines are the the 2x2a error-snakes. High-extinction 
SNe (blue crosses) are above the maximum possible intrinsic colors 
calculated from the SNe in passive galaxies and were not included 
in the fits. CC-SNe (red triangles) are also shown. 



2.2. Luminosity- color- stretch 

We next derive a relation between the absolute peak 
magnitude as a function of stretch and color. To limit 
the effects of peculiar velocities, we require SNe in our 
training sample to have either zcmb > 0.01 (CMB frame 
redshift) or a known distance to the host taken from the 
NED database. 

We fit simultaneously magnitude, stretch and color 
to different models. The best one is again bilinear in 
both, stretch and color. The absolute magnitude shows 
a steeper slope as a function of color for redder objects 
changing at c ~ 0.07 (upper panel, Figure [3]). Against 
stretch, the magnitude values (lower panel, Figure[3]) are 
also well fit by a bilinear model with a shallower slope 
for the low-s stretch range. The effects in the differing 
slopes are not as strong as for the color-stretch relation 
but they will make the posterior analysis more precise. 
The errors in the absolute magnitude (from the appar- 
ent magnitude and redshift errors) and in the stretch and 
color are included in the fit. The magnitude-stretch-color 
relation can be summarized as follows: 



match 



M B = M + a [lj3] (s - 1) + (3 [aM c + M^jf 



(2) 
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withM = -19.11 ±0.01, and where a x = -1.19 ±0.03 
and «2 = —0.54 ± 0.09 are respectively used for stretch 
values greater and lower than the transition value of s = 
0.82 ± 0.02. f3 a = 3.00 ± 0.02 and (3 b = 2.15 ± 0.05 are 
used for values of c greater and lower than 0.07 ± 0.01, 
respectively. M^ atch = -0.18(0!! - a 2 ) and M b " latch = 
0.07(/3 a — jib) unite the fit at the changing values. This 
model results in a slightly better fit (xf, — 19.6) than 
other models (e. g. exponential -in stretch and color- 
with xl = 19.9). 
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Fig. 3. — Absolute magnitude (corrected for stretch from Eq. [2] 
i. e. MJf corr = M B - £*ri, 2 ] (s- 1) - M™ atch ) vs color (upper) and 
absolute magnitude (corrected for color, i. e. M B ~ corr = — 
/^[a,b] c — M™ atch ) vs stretch (Zoiuer) for the training sample. The 
solid lines represent the best fit model, a two-piece linear model for 
each case. The dashed lines are respectively the 2 X 2a error-snakes. 
CC-SNe (red triangles) are also shown. 

We include an asymmetric error-snake on these rela- 
tions, shown by the dotted lines in Figures [2] and 
These delimit (twice) the region within which 95% of 
all SNe la in each bin of 0.1 in color and 0.07 in stretch 
are found around the best model. The width of the re- 



gion has been doubled to allow for the larger photomet- 
ric scatter of the high-z SNe, which will be important 
later in £|3.2[ This region also covers any uncertainties in 
the relations modeled in this section and allows for scat- 
ter in the color of the SNe. Requiring all objects to be 
inside those error-snakes, i.e. color-stretch, magnitude- 
color and magnitude-stretch, we find that ~ 90% of the 
normal and low-s SNe la pass the cut (see Table[3]). How- 
ever, only 2 contaminants (all SNe Ib/c) survive the same 
cull, a ~ 6% of the total of contaminants. We find this 
tradeoff to be beneficial. Reducing the contamination is 
useful when limiting the systematic errors as discussed 
in 



3. HIGH-Z LOW-STRETCH CANDIDATE SELECTION 

The search component of the SNLS, carried out using 
the Canada- France-Hawaii Telescope (CFHT), is a wide 
field imaging survey well-suited to the detection of SNe l a 
in the redshift range 0.1 < z < 1.0 (lAstier et al.ll2006f ). 
It is a rolling search where four 1° x 1° fields (Dl to D4) 
are imaged every few days, so that a well sampled LC 
is obtained in four filters (g m r M ^mzm) for each discov- 
ered supernova. Good candidates are followed-up spec- 
troscopically. Because they are faint, less frequent, and 
not vital for cosmological studies, typical subluminous 
SNe la are not pr ioritize d in the spectroscop ic follow- 
up (|Sullivan et aij|2006at iPerrett et al.ll2010bl) . We can 
still use the well-sampled SNLS multi- wavelength LCs to 
photometrically identify these objects. The method used 
resembles the one fo r the normal-stretch population of 
IPerrett et all ()2010aL hereafter P10), although with ap- 
propriate adjustments for the low-s objects. This makes 
a direct comparison with the different stretch popula- 
tions possible. 

During the first 4 years of operation, up to July 2007, 
the candidate database consists of more than 5000 photo- 
metric transients including AGN, non-SNe and CC-SNe. 
As in P 10, candidates in masked areas of the deep stacks 
are removed, as well as objects that lack the necessary LC 
coverage for a good photometric identification and pa- 
rameterization. These observational criteria reduce the 
sample to half and are the same ones used for the SNLS 
training sample in |J2] at least one §m observation be- 
tween — 15 and 5 days of maximum, at least one Tm and 
one im observations between —15 and ±2.5 days of max- 
imum, at least one r M and one im observations between 
—9 and ±7 days of maximum, and at least one tm or one 
%m observation between +5 and ±20 days of maximum. 

3.1. Redshift estimation and LC fit 

In P10, accurate photometric redshifts from the SNe 
were obtained with the use of two normal SN la LC fit- 
ters. Their approach requires a good spectroscopic sam- 
ple for testing the method and the precision of the final 
photometric redshifts. We do not possess any 91bg-like 
SN la spectrum at high redshift with clear Ti II fea- 
tures, only 3 low-s objects at the border of our defintion, 
s < 0.8, that have spectra of rather normal SNe la. The 
high degree of degeneracy in the LC fits, particularly 
with stretch and redshift, which compete to make an ob- 
ject appear faint, demand an unbiased redshift estimate 
that will lead to appropriate LC parameters like stretch 
and color. 
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Several hundred objects of the 4-year dataset possess 
spectroscopic redshifts (spec-z) from either the SN or 
from the host galaxy. For the SN-derived redshifts, the 
SN type has bee n determined from t he spectra in the ma- 
jority of cases (IHowell et all 12003 IB render et all 12001: 
lEllis et aLlfeOOaiBalland et al.ll2009l ). Most of these are 
normal SNe la (413), although a good sample of CC 
objects is also available (78 ). The host spec- z include 
data from DEEP/D EEP2 (iDavis et al l l200l. V VDS 
(ILe Fevre et al.ll2005lV zCOS MOS (ILil lv et al.ll2007r i and 
VIMOS (jJonsson et al1l2010( ). Many candidates without 
spectroscopic redshifts have host photometric redshifts 
from the SNLS broad-band photometry either using the 
spectr al connectivity analysis method of Freeman et al. 
( 2009), or, when no match is fou nd, through fit s to spec - 
tral energy distributions as in ISullivan et afl ([2006b). 
The resulting host photo-z, Zhostphot, are very good to 
z ~ 0.75 with ' Az| , = 0.042 and a low (~ 5%) catas- 
trophic failure rate ( Jjas^oatphotj > 9.15). 

Fortunately, some of those SNe la with catastrophic 
redshifts can be identified with the help of the two 
LC fitters of P10 adj usted for low-s SN: estimate_sn 
(jSullivan et al.ll2006ali and the aforementioned modified 
version of SiFTO (see ^j. Their use in the multi-step 
process explained below provides good quality-of-fit pa- 
rameters capable of distinguishing SNe la from CC-SNe, 
but also offer an independent mechanism for obtaining 
photometric redshifts using the SN, zsNphot, in cases 
where host redshifts are lacking. 

Once the redshift is known, SiFTO is used to obtain 
the stretch and color of each SN. Tests on the local sam- 
ple (Sj2| show that it works well without imposing any 
stretch-color model. As shown in Figure 2 of P10, the 
derived stretch is largely independent of the redshift be- 
cause the rise and fall of SNe la are faster at shorter 
wavelengths, which largely compensates for redshift un- 
certainties. In Figure [U we demonstrate that this behav- 
ior extends to low-s objects. Shifting the redshift -either 
photometric or spectroscopic- for each SN in our final 
low-s sample (see §3.2[) by as much as Az = ±0.3, we 
obtain a stretch variation lower than ~ s ± 0.1. While 
this shows more dependence on redshift than the normal- 
stretch population, it still allows us to obtain accurate 
stretch estimates even in the presence of moderate red- 
shift uncertainties. 
The estimate_ sn rout ine is described in more detail 
(p006ah . 



in ISullivan et al.1 ()2006al ) . It performs a grid search on 
the x 2 01 a nt t° SN photometry in flux space to a five- 
parameter model, including redshift. For this analysis, 
we modified it to include low-s objects using the color- 
stretch relations of Eq. [TJ and the magnitude relation of 
Eq. [2] This fit depends on the assumed values of the 
cosmological parameters. 

We calculated the SN photo-z for all objects in a simi- 
lar fashion to P10 using the following multi-step process: 
First, we use estimate_sn to get a preliminary redshift 
estimate. We then use this redshift as a fixed input to a 
SiFTO fit to obtain an improved stretch measurement. 
Finally, we feed back this stretch to a final estimate_sn 
fit as a fixed parameter to obtain a more accurate zsNphot- 
We use this SN photo-z for SNe that do not have any 
other redshift (spectroscopic or host photometric). 




Fig. 4. — Upper: The change in fit stretch due to the deliberate 
variation in redshift for the low-s final sample of i]3.2l The blue 
squares represent the spectroscopically-confirmed low-s SNe la. 
The solid black lines are the standard deviation above the me- 
dian As. Lower: The average output stretch for each SN la as a 
function of its stretch (at zero redshift offset). The error bars are 
the change due to the variation in redshift for each SN la. 

Objects with Zhostphot that have "bad fits" , i.e. that do 
not pass all cuts of £13.21 are not discarded, but instead 
the cuts are re-applied using zsNphot- If this also fails, 
then the object is rejected. This offers some protection 
against catastrophic host photo-z failures. This recov- 
ers 29 confirmed SNe la. All known contaminants fail 
with both photo-z values, so we are confident that this 
procedure still rejects non-las. 

3.2. Final candidates 

We perform final SiFTO (and low-s-modificd 
estimate_sn) fits to each object to obtain the pa- 
rameters of each SN and quality-of-fit values. A 
series of cuts are then applied to screen out objects 
with poorly estimated parameters or those unlikely 
to be SNe la. A summary of the applied culls and 
the number of objects surviving each of them can be 
found in Table SJ The \ 2 culls are similar to those of 
P10, although we weaken the requirements slightly for 
s < 0.8 objects based on the study of the low-z sample: 
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reduced x 2 s 01 Xe S timate_sn < 13 for the overall fit, and 
for individual passbands Xc S timatcjsn( r A;f ) < H an< i 
Xe S timat OJ3 „(u/) < 16, and a SiFTO cull of xIfto < 11 • 
Additionally, we require the objects to be inside the 
error-snakes trained in $2] with low-z low-s SNe in order 
to constrain their parameters (Figure EJ}. These cut 
values are derived empirically by trying to minimize the 
number of known contaminants that pass and minimize 
the number of rejected confirmed SNe la. 




*-spec 

Fig. 5. — Photometric redshifts obtained from the host or from 
the LC fits compared to spectroscopic redshifts from the SN (open 
circles) or the host (filled circles). Purple stars represent s < 0.8 
objects. The solid line shows 1:1 agreement and dashed lines a 10% 
uncertainty in (1 + z Bpcc ). The photo-z errors are either from the 
host or from estimate_sn. 



3.2.1. At z < 0.6 

Up to z = 0.6, the final sample contains 161 objects, 
of which 126 have spectroscopic redshifts. They span the 
stretch range 0.6 < s < 1.2, of which 18 have s < 0.8. 
Six of these have known spectroscopic redshifts: three 
are confirmed SNe la with s ~ 0.8 with no outstanding 
91bg-like feature in their spectra, and the other three 
have spectroscopic redshifts but have too low of a con- 
trast with the underlying host galaxy to conclusively es- 
tablish the spectroscopic SN type. The final parame- 
ter distributions are shown in Figure [7] and Figure [8l 
The low-s photometric sample has median stretch and 
color of (s) = 0.74 and (c) = 0.24 (as compared to 
(s) = 1.00 and (c) = 0.00 for the normal-s population). 
The mean (Amag) = 0.76 ± 0.69 clearly shows that our 
low-s sample is not as extreme as the characteristic 91bg 
of Amag ~ 2mag. 

The complete final set of photo-z contains 82% with 
host photo-z and 18% with SN photo-z. Figure [5] shows 
the comparison of photometric and spectroscopic red- 
shifts for all SN la candidates with spectroscopic redshift 
either from the SN or the host. The median precision is 

_[ Az — 7 = 0.030. The photo-z are more accurate to 

(, -1 \Z spec ) 



z ~ 0.6 (|Az|/(l + z spec ) — 0.025), where our rate study 
will concentrate. 

The final errors in photo-z are obtained from Figure [5] 
by taking the standard deviation around the median as 
a function of z. The stretch and color errors are pro- 
vided by SiFTO, which requires a known redshift input. 
Stretch and color variations due to redshift uncertainties 
will be taken into account with a Monte Carlo simulation 
in §[H2j 

3.2.2. At z> 0.6 

The detection incompleteness at 0.6 < z < 1.0 in- 
creases rapidly for low-s SNe la, as will be shown in next 
section. We still find a significant population of 44 can- 
didates (11 with spec-z from the host and 5 with low 
signal-to-noise spec-z from the SN, 2 of which are con- 
firmed SNe la of s ~ 0.8 and no 91bg features) with 
median stretch, (s) = 0.76, and color, (c) = —0.02. The 
average color is extremely blue compared to the median 
color at z < 0.6 of 0.24 confirming that we are only ob- 
serving the blue tail of the low-s distribution. This can 
be seen in Figure^ They are not very faint objects, even 
for low-stretch, as shown in Figure [5] The lack of char- 
acteristic low-s red objects makes it difficult to reliably 
estimate the low-s rate at z > 0.6 and proper corrections 
need to be made (see § 14. 2[) . 

Additionally, we do find a presence of blue low-s can- 
didates at high-z which are absent at lower-z. Such blue 
objects with c < —0.05 are neither found in the nearby 
nor in the z < 0.6 SNLS samples, but represent ~ 30% 
(14) of all z > 0.6 low-s objects. We do not expect all of 
these objects to be CC-SNe, which could potentially con- 
taminate the low-s sample with their characteristic bluer 
colors. As will be shown later, besides going through our 
SN la pipeline, the environment of these candidates is 
not characteristic of CC-SNe. 

A second possibility is that this group represents an- 
other class of tr ansients that are fa st and blue. The 
object found by iPerets et al.1 ((2010) is faint and fast- 
evolving but with clear different spectroscopic features 
and bluer colors than 91bg-like objects. It is hosted in an 
old stellar population as well as most of our candidates. 
However, our objects are not as faint (Figure Et- On 
the ot her h and, the transients fou nd by |Poziiaiisk i"et al.l 
(|2010f) and iKasliwal et al.l (|2010j ), suggested SN .la-like 
helium detonations of a WD, are blue and of similar 
brightness than ours, but evolve faster. When fit with 
our templates, we obtain very poor fits with s = 0.44 
and s = 0.33, respectively, lower than our lower stretch 
limit for subluminous SNe la. 

A tempting explanation is that the input photometric 
redshift for the LC fit is consistently wrong and affects 
the obtained color and stretch. This hypothesis is ruled 
out because 7 of these objects have spectroscopic redshift 
and the rest (except one) have photometric redshift from 
the host (with mean photo-z error of ~ 0.08), which do 
not show any significant systematic. We have seen the 
small effect in the measured stretch when the redshift is 
varied (Figure 01 ■ We now show the effect in the mag- 
nitude and color of photometric candidates in Figure [TU1 
Their variation due to changes in photo-z cannot account 
for all the bluer colors but for some of them -only if the 
redshifts were all consistently overestimated by +0.05, 
the colors would be redder by +0.08 (and magnitudes 
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Color-stretch (left), magnitude-color (middle) and magnitude-stretch (right) for all SNLS objects passing the observation 
criteria (open points) and the final SN la photometric z < 0.6 sample after all cuts of Table [4] (filled points). Spectroscopically-confirmed 
SNe la (blue squares) and non-la (red triangles) are shown. The lines are the best-fit relations and error-snakes found in [[2] for the local 
training sample, and served as additional cuts. 
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Fig. 7. — Final distributions in redshift (left), stretch (center) and color (right) of the z < 0.6 SN la photometric sample including low-s. 
The white histogram is the whole population, the gray is the spectroscopically-confirmed SNe la and the black the low-s candidates. 
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Fig. 8. — Distribution in Amag of the final low-s candidate sam- 
ple. Amag = M + 19. If is the difference between the absolute 
magnitude of each SN compared to the magnitude — f9.flmag of 
a "standard" SN la at s = 1 and c = according to eq. [2] The 
shaded histogram shows only the low-s SNe la at z < 0.6. 



fainter by +0.17, whereas stretches would only change 
by -0.004). 



Fig. 9. — Color as a function of redshift for the final photometric 
SN la sample. Normal-s (black circles) and low-s (purple stars) 
are highlighted. 



Photometric errors in stretch and color get larger as we 
probe higher redshifts. Most of the candidates found at 
these redshifts have stretch close to the limiting defini- 
tion of s = 0.8, so that they are consistent with being fast 
normal-s SNe la. Only 3 candidates have s + s crr or < 0.8, 
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Fig. 10. — Absolute magnitude vs color for the low-s SN la sample 
with photo-z at 2 < 0.6 (big open pink stars) and at z > 0.6 
(small filled purple stars). The dots show the normal- s sample. 
The upward and downward arrows show the new magnitudes and 
colors when varying the photo-2 by +0.05 and —0.05 respectively. 
The lines unite the same object at the three redshifts. The lower 
right error bars represent the typical (mean) photometric errors. 

one of them with s < 0.6. Of these, two also have 
c + Cerror < —0.05, inconsistent with typical low-s col- 
ors, and will be further discussed in 14.31 The rest can 
therefore be explained as the blue tail of s ~ 0.8 objects 
while the redder counterpart remains unseen due to se- 
lection effects. In the remainder of this work, we treat 
this group as such and correct for this selection bias when 
calculating the rates. 

4. LOW-STRETCH SN IA RATE EVOLUTION 

With a set of low-s candidates, we proceed to calcu- 
late a volumetric rate as a function of redshift. As in 
P10, the volumetric rate, rb(zb), for a given redshift bin, 
represented here by the middle of the redshift bin, Zb, is 
given by the sum over all Nb candidates inside the bin: 

1 (1 + Zj) 

2—1 

where e is the individual SN detection efficiency de- 
pending on the characteristics of each object (like red- 
shift, stretch and color) and AT is the sampling time ac- 
cessible to find the object, divided by (1 + z) to account 
for time dilation. V b is the comoving volume contained 
in the redshift range of the bin. 

4.1. Detection efficiencies 

The detection efficiencies ti account for observing bi- 
ases measuring the rate. They depend on numerous fac- 
tors and are specified as a function of the SN charac- 
teristics such as redshift, stretch and color (along with 
field and year of observat ion). The effic i encies used here 
are described in P10 and lPerrett et al.l (|2010b| ). and are 
based on a Monte-Carlo study of 2.4 million artificial 
SNe la injected into the actual SNLS images. This ap- 
proach avoids the need to make assumptions about tem- 
poral and spatial data coverage. 

The artifical SNe are generated using uniform distribu- 
tions of redshift and stretch (down to s = 0.5), and the 



color is determined using a relation like Eq. [TJ including 
an intrinsic scatter around the stretch-color relation (a 
Gaussian noise with dispersion a — 0.04 added indepen- 
dently to U and V) . The peak magnitudes are generated 
from the luminosity-stretch-color relation, as in Eq. [2j 
and the redshift and assumed cosmology, plus an addi- 
tional unmodeled magnitude dispersion from a Gaussian 
distribution with an intrinsic scatter of (T mf = 0.15. The 
LC is generated from integrating, at a given epoch and 
through the filter functions, the SN la spectral templates 
scaled for the B-magnitude and corrected for U — B and 
B — V colors, Milky Way and host extinction. The ar- 
tificial SN is then inserted into each iM image obtained 
by the SNLS. The images are then processed through 
the same pipeline used in the real-time SN discovery, 
and the recovered fraction is used to derive the efficiency 
as a function of various parameters, including redshift, 
stretch, color, observi ng season and field (see a more 
detailed description in iPerrett et al.l l2010bD . The low- 
s candidates, which are faint, red and have a rapid rise 
and fall, are more difficult to detect and therefore have a 
lower efficiency, so are accordingly given a higher weight 
when calculating the rates. 

The efficiencies are calculated in bins of width Az — 
As = Ac = 0.1. The measured efficiencies, linearly in- 
terpolated according to each SN, are shown in Figure [TT1 
It can be seen that the mean efficiencies are lower for 
the low-s sample up to z ~ 0.6. Above z — 0.6, the 
mean efficiencies seem to approach the normal- s sample. 
At these redshifts they become low for all but the bluest 
low-s objects. This leads to a higher fraction of low-s 
blue objects with very few red ones recovered, as found 
in § 13.21 We show the expected efficiency distribution 
for the missing red, low-s SNe la of an unevolving color 
distribution as green squares in Figure [TTJ based on the 
observed low-s color distribution at low redshift from iJTl 
The rates at z > 0.6 will need to be corrected for these 
missing objects (see next section). We note that the effi- 
ciencies shown in Figure [TT] are the raw efficiencies -the 
recovery fraction from the simulated SNe- without the 
observational constraints of section § [3] nor the sampling 
time correction that depends on each field and season. 
These corrections can make the actual efficiencies much 
lower (see P10). For instance, one of our low-s candidates 
at high-z with z = 0.80, s = 0.63 and c = 0.24 has a raw 
efficiency of 0.54 but with observational constraints and 
sampling time this efficiency is only 0.18, which explains 
why these objects are hard to observe. 

4.2. SN la Rates and Monte-Carlo simulations 

We calculate a rate evolution, r meas , shown in Column 
5 of Table[5]solely based on Eq.[3]and with the efficiencies 
from § 14.11 We use bins of Az = 0.2 due to the low sam- 
ple statistics (the last z-bin is only 0.1 wide as no low-s 
candidate exists beyond z = 1). The measured errors 
in stretch and color from the LC fit, and especially the 
errors in the photometric redshift of each candidate, can 
lead to uncertainties in the rate measurement. We ac- 
count for these via a Monte-Carlo (MC) simulation that 
randomly shifts the photometric redshift according to its 
error (for those SNe la with photo-2). For each redshift, 
a corresponding LC fit provides new parameters and its 
covariance matrix. This latter is used to generate a new 
random set of stretch and color using standard tech- 
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Fig. 11. — Efficiencies calculated by P10 for the SN la photomet- 
ric sample for normal-s (black circles) and low-s (purple stars) pop- 
ulations. For objects where no direct efficiency was found (small 
filled stars), the closest match in color was used. Mean values 
in redshift bins for normal-s and low-s populations (big filled cir- 
cles and stars) are shown. Additionally, simulated low-s objects for 
2 > 0.6 bins (green squares) based on the z < 0.6 color distribution 
and respective mean (green dashed line) are shown for comparison. 
Note that these efficiencies do not include the observation criteria 
of section § \3\ and sampling time correction of the different fields 
observed. 



niques (| James fc RooslfT975l) . We recalculate the rates 
from Eq. [3] for each iteration based on the randomized 
values, and use the median value over the 1000 realiza- 
tions to determine the final rate, tmc- We estimate the 
error from the variance of the rates in each bin (and call 
them "MC" errors shown in Column 6 of Table [5] and 
Figure [T2|) and add it in quadrature to the weighted la 
Poisson statistical errors in Figure [12] (confidence levels 
calculated using the Clop per-Pearson method used for 
example in iGehrelsl 119861 ) . The error simulation differs 
from P10, as the SNLS does not possess a known spectro- 
scopic sample of low-s SNe la from which we could obtain 
parameter distributions. Instead, we treat each SN indi- 
vidually by randomly drawing parameters according to 
its errors and covariances. Although the precise details 
of our simulation differ from P10, a SN la rate evolution 
within the errors is obtained. 

Additionally, we have seen in the previous section that 
we do not find red objects in the higher- z bins. This 
selection bias results in underestimated rate values at 
high-z, which are not corrected by the efficiencies of the 
bluer objects found. There is a lack of objects in certain 
red color bins, so that no efficiency can account for them. 
To measure this effect, we take the low-s distribution of 
color in an underestimated high-z bin (for each MC iter- 
ation), we normalize this distribution to the known color 
distribution at z < 0.6 (right Figure [7]), which has al- 
ready been corrected for the efficiencies of each SN, and 
then generate fake SNe in the missing color bins. The 
normalization is done with the two bluest matching bins. 
Figure IT51 shows an example for 0.7 < z < 0.9. The green 
open histogram is the final corrected histogram from the 
initial observed diagonal-shaded distribution. The cor- 



recting factor for this case is /fc^orr = 1 + jVb ^ kc ^ = 1.9, 
where A";, (fake) is the number of generated SNe in the bin 
and Ni, is the observed number, so that the corrected rate 
IS V}} corr = n X /b,corr- The final correcting factors from 
the median of all MC iterationns, summed Poisson (from 
the number of objects used in the bins used to normalize 
both distributions) and MC errors, as well as the result- 
ing corrected rates are shown in the last two columns of 
Table [5] In Figure [T^l the corrected rates are shown as 
stars. 

Finally, errors in the efficiencies can affect the rate val- 
ues. As explained in P10, those come mainly from sys- 
tematic errors in the underlying assumed Amag distri- 
bution of SNe la, the distribution of residuals between 
the absolute magnitude of a fiducial s = 1 SN la and 
the magnitudes of observed SNe la after correction for 
stretch and color. To check for this systematic error, wc 
allow for random efficiencies in the MC simulation cal- 
culated from a range of dint = 0.12 — 0.15 in the Amag 
distribution, obtaining results that agree very well with 
the rates of Figure [T3] 

4.3. Contamination uncertainties 

An important systematic error comes from core- 
collapse SNe mimicking the LC of low-s SNe la. In P10, a 
missclassification fraction of less than < 3% was found for 
the normal-s population. We expect this fraction to be 
higher for low-s SNe la due to two effects. First, they are 
dimmer, and more easily confused with the intrinsically 
fainter CC-SNe at a given redshift. Second, as we probe 
higher redshift ranges, CC-SNe at low-z might masquer- 
ade as low-s SNe la at higher-z since we do not have 
spectroscopic redshifts for all candidates. SNe II have 
characteristic LCs that are very different from SNe la 
and present less danger of contamination. SNe Ib/c, on 
the other hand, re present a smaller f r action of CC-SNe, ~ 
30% a ccording to ICappellaro et ail (119991); iSmartt et alj 
(|2009ft and 25% according to |Lj_etaL| (J2010aj), but are 
more similar to SNe la and will be the most probable 
contaminants. 

The first error -where the photo- z is correct, yet the 
fit of a SN la template manages to pass the quality- 
of-fit cuts- can be examined using the training sam- 
ple (as in Ej2j) . We found a 6% CC-SN missclassifica- 
tion rate, all SNe Ib/c. This fraction can be propa- 
gated to higher redshift using the SN la rate evolution 
of P10, Rj a {z) cx (1 + z) ata with a Ia ~ 1.91, and a 
CC-S N rate evolution reflecting the star formation rate 
from lHopkins &: Beacoml (]2006f ). Rcc (l + z) acc with 
acc — 3.6. We assum e a ratio of SNe la to CC-SNe of 
~ 0.22 at z = 0.3 from iBazin et al.l (|2009f ) . which trans- 
lates to a local ratio of ~ 0.39 - as compared t o ~ 0.48 
found b v lCappellaro et al.1 (JT999) and 0.24 from lLi et"afl 
(|2010al) -. The percentage fraction of contamination as 
a function of z, fcont(z), can then be described as: 



/cont(z) 



passcont 



A pass (z) l + 0.39r(l + z) 



ai a — QCC 



(4) 

where r = 0.88/0.06 is the cull passing efficiency (the 
probability of a SN passing all selection criteria) ratio 
of la with respect to CC, assuming it does not evolve 
with z. This equation allows us to get a crude estimate 
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Fig. 12. — Volumetric rate of low-s SNe la as a function of redshift (green open points) compared to the normal-s rate (blue squared 
points) of P10. The total summed rate (blue filled points) is shown with a little z-offset for clarity. The stars represent the corrected high-z 
low-s rate calculated from real and simulated objects drawn from a distribution scaled to the z < 0.6 color distribution. The filled triangle 
is the total nearby SN la rate from Li ct al. (2010b) and the open triangle is the low-s rate estimated from a s < 0.8 fraction of it. The 
solid lines are the best A + B model and uncertainties found in P10. The dashed lines are the best A-model and uncertainties found in this 
study through a fit of only SNLS z < 0.6 bins: A = (1.05 ± 0.28) X 10" 14 SNeyr-iM" 1 . 



of the number of contaminants that would pass all cuts, 
iVpasscont(z)- By dividing the expected fraction of con- 
taminants per z-bin by the mean weight of low-s SNe 
found therein (per MC simulation), we translate this into 
a rate error lower than 20% for z < 0.6 up to as high as 
30% at z ~ 0.9. The estimates are sensitive to the power 
law of the rate evolution, the initial fraction of SNe la to 
CC-SNe and the local missclassification fraction of 6%. 
Increasing this latter fraction to 10 % or usi n g SN l a to 
CC-SN local ratio of 0.241+^ by ILi et all (I2010al) in- 
stead, we obtain contamination percentages up to 30% 
for z < 0.6 and 40% at z ~ 0.9. This latter rate error 
from contamination is added in Table [5] and we believe it 
to be an upper limit. The contamination error is smaller 
than the combined statistical and MC errors at z < 0.6 
but becomes comparable at higher-z. One should note 
that the missclassification error increases with redshift 
and would only make the rate higher than the actual 



The second source of contamination comes from uncer- 
tainties in the photo-z. Incorrect redshift estimates can 
cause low-z CC-SNe to simulate low-s SNe la at high- 
z. This effect can be tested by running all spectroscopic 
objects through our selection criteria with their calcu- 
lated photo-z instead of the spec-z. By doing this, we 
retain most of our candidates with a similar fit stretch 
and, most importantly, we still reject all the CC-SNe as 
contaminants (even though their photo-z was calculated 
with SN la templates and can occasionally be very dif- 
ferent from the spec-z). 

Perhaps the most concerning signature of contamina- 
tion is the low-s SN la sample at z > 0.6. These have 
rather blue colors and could be CC-SNe instead. As an 
additional test, we analyzed the host galaxy properties 
(Figure 03]). In general these objects are found in mas- 
sive and passive galaxies, not characteristic of CC-SNe. 
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Fig. 13. — Example of how the observed color distribution for 
low-s at 0.7 < z < 0.9 (diagonal-shaded) is corrected to match the 
efficiency-corrected low-s color distribution at z < 0.6 (shaded). 
The final corrected color distribution for that z-bin (open his- 
togram) has /(, = 1.88 more SNe than the original observed one. 

This provides evidence that most of these candidates, in- 
cluding the z > 0.6 population, are likely the actual blue 
tail of the low-s SN la distribution, as mentioned earlier. 

Nevertheless, two of these objects, 05D3hp (with s = 
0.59 ± 0.04 and c = -0.16 ± 0.09) and 07D2bc (with 
s = 0.72 ± 0.04 and c = -0.13 ± 0.03), have unusually 
blue colors at very low-stretch. The host of 05D3hp is 
small (M = 10 7 ' 9 M Q ) and star-forming (SFR/mass = 
10 _9 8 A/oyr _1 ) and points rather towards a CC event 
imitating a low-s SN la. 07D2bc happened also in a 
star-forming (SFR/mass = 10~ 100 -MQyr _1 ) but rather 
massive galaxy (M = lO n -°M ). If these two objects 
are CC-events, they would constitute a ~ 5% contam- 
ination fraction. Although it is dangerous to assume 
the local star formation from general galaxy properties, 
if we were to assume that all events in "burst" hosts 
with SFR/mass < — 9.5M Q yr _1 are contaminants, then 
~ 30% of the z > 0.6 objects are missclassified (~ 25% 
at z < 0.6), in agreement with our generous estimates. 

5. DISCUSSION 

We now turn to the interpretation of the low-s SN la 
rate in the light of suggested rate models, and com- 
pare it to different stretch populations. We investi- 
ga te the empirical t wo-c ompone nt rate evolution model 
of IMannucci et all (|2005| ) and iScannapieco &: Bildstenl 
( 2005f) but also study simple delay-time dis tributions 
( Dahlen fc Franssonl[l999l : IStroker et al.ll200"l for SN la 
progenitors. A local low-s rate estimate based on nearby 
studies helps constrain the models. 

5.1. Low-z subluminous rate 
Recently, iLi et al.l ([2010b[ ) have calculated a low- 
redshift SN la rate and a corresponding 17.9lg 2% lrac " 
tion of 91bg-like objects f)Li et al.l [20 10a). To obtain a 
rough local rate estimate of low-s SNe la defined accord- 
ing to our study, i.e. s < 0.8, we use thei r luminosity- 
function sample (Table 3 of ILi et aLll2010al ) and directly 
fit objects for which we have photometric data (48 of 
74) with SiFTO. For the rest we fit instead their clos- 
est i?-band LC template match (column 11). Using the 



completeness-corrected numbers (column 9), we obtain 
16.4 objects with s < 0.8, or a 21.9^1% fraction. This 
fraction is complete within a volume of 80 Mpc and trans- 
lates to a low-s rate of (6 .3 ± 1.8) x 10~ 6 y r- 1 Mpc~ 3 ^ n 
based on the local rate of ILi et all (|2010q ) (adjusted for 
^70 = 1) an d including the uncertainty from the Poisson 
error. We note that by using the whole low-z training 
sample of Sj2j the 41 low-s objects correspond to a 22% 
fraction of the total, in good agreement. 

This rate estimate is based on small sample statistics 
and suffers from the biases that arise from a distance- 
limited host-targeted survey, where the rate of SNe is 
converted to a volumetric rate via the galaxy luminosity 
function. This can lead to systematic errors, as brighter 
galaxies may contain more SNe la of a certain type (as is 
expected for low-s). However, a local low-s rate estimate 
is a valuable tool for comparison to the high-z rates. 

5.2. Power-law evolution 

The low-s SN la rate, as opposed to the normal-s pop- 
ulation, does not seem to increase with redshift and fol- 
low the star formation history (SFH) trend found in P10 
(Figure IT2|) . If we fit a simple power-law, r oc (1 + z) a , 
to the SNLS rate to z = 0.6 and the low-z estimate, 
we obtain a — —0.90 ± 1.38, smaller than the value 
of the normal-s population found by P10 (a n0 rmai = 
1.91 ±0.24). This means that an evolution similar to the 
normal-s population is rejected with ~ 94% probability. 
If we include the corrected z > 0.6 SNLS points, we find a 
higher slope more consistent with zero: a = — 0.34±1.41. 

5.3. Delayed component model 

It is known that SNe la occur in both star form- 
ing and non-star forming regions, whereas low-s objects 
are usually found in passive ones. In the model o f 
IMannucci et al.l (|2005l) ; IScannapieco fc Bildstenl (|2005fl , 
two populations are considered: a "prompt" component 
for SNe exploding in 0.1 Gyr or less that is proportional 
to the star formation history, M(z) and a "delayed" com- 
ponent for events with higher delay times that is propor- 
tional to the host mass, M(z). It is described as: 

SNR I& (z) = AM{z) + BM{z), (5) 

where A and B are scale factors for both components. 
The delayed model for low-s SNe la is supported when 
one looks at the host properties of the low-s sample. 
CFHT Legacy Survey host galaxy broadband photom- 
etry can be fitted to s everal template galaxy SEDs us- 
ing P EGASE/ZPEG dLe Borgne fc Rocca-Volmerangd 
2002) to determine average galaxy star formation 
rates (SFRs), stellar masses and photometric redshifts 
([Sullivan et all l2006b[ ). The properties of SNe la are 
compared to their host type (see Figure fl4j) to inves- 
tigate the behavior of low-s objects in contrast to the 
normal-s population. 

We use the specific star formation rate (sSFR), the 
star formation per unit stellar mass, as the indicator 
of star formation activity, and separate the host galaxy 
sample in passive galaxies, log(sSFR(yr -1 )) < —11.5, 
weak star formers, —11.5 < log(sSFR(yr -1 )) < —9.5, 
and strong star formers log(sSFR(yr -1 )) > —9.5, where 
SFR refers to the average star formation over the last 
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Fig. 14. — Host properties for z < 0.6 (upper) and z > 0.6 (lower) SN la candidates. Left: Normalized specific star formation (star 
formation per unit mass) distribution divided in normal-s (blue) and low-s (red) populations. The probability that the two populations 
are drawn from the same distribution (KS-statistic) is less than ~ 0.5% for both redshift ranges. Right: Normalized host stellar mass 
distribution per stretch range, normal-s (blue), and low-s (red). The KS-statistic is below ~ 7% for both redshift ranges. 



0.5 Gyr. This approach ignores important internal vari- 
ations of SF in the galaxy, as shown in Ras kin et al.1 
(20093), and provides only a general estimate of the SF. 
The sSFR distribution for two different stretch ranges 
(left Figure [14]) , normal-s and low-s, provides clear ev- 
idence that the low-s populatio n is less dependent on 
star formation. iNeill et all (2009) find that no local low- 
s SN la is found in a host with log(sSFR(yr^ 1 )) > -9.5. 
Here, we corroborate their results for high-z objects find- 
ing that only 3 (17%) at z < 0.6, and 11 (25%) at higher 
redshift, have higher mass- weighted star formation rates. 

Moreover, the host mass distribution for these stretch 
ranges (Figure [Til right) shows clearly different distribu- 
tions, in which low-s SNe la occur almost solely in m as- 
sive hosts. We extend the results of lNeill et al.l (|2009l ) to 
the distant universe: low-s SNe la are found mostly in 
hosts with masses above ~ 10 9 ' 5 M Pi . For all cases, the 



distribution of host-galaxy mass for the two populations, 
normal-s and low-s, have a low two-sided Kolmogorov- 
Smirnov (KS) statistic, indicating a small probability 
that they are randomly drawn from the same distribu- 
tion. 

The above relations of SN la properties with host char- 
acteristics demonstrate that previous findings extend to 
the extreme of the SN la population at z > 0.1, and at 
the same time are a confirmation that the observed pho- 
tometric low-s sample obtained is real and not strongly 
contaminated from core-collapse events that are only 
found in late-type galaxies. 

They also provide additional support to model the 
low-s SN la rate as just an A-component. Fitting the 
latter to the local and SNLS and rates up to z = 0.6 
provides A = (1.05 ± 0.28) x 10" 14 SNeyr^M" 1 and 
is shown in Figure 1121 If we fit both components, we 
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obtain A = (1.29 ± 0.66) x 10~ 14 SNeyr-^g 1 and 
B = (-0.27 ±0.64) x lO^SNeyr-^Mgyr- 1 )" 1 . When 
we include the corrected SNLS z > 0.6 points, we ob- 
tain A = (1.14 ± 0.61) x lO-^SNeyr-iM" 1 and B = 
(0.03 ±0.62) x lO^SNeyr-^Meyr- 1 )- 1 . In the latter 
two cases, B is close to zero (B < is an unphysical 
result) suggesting that the mass component dominates 
over the component proportional to the star formation 
history. 

Comparing our delayed component (for z < 0.6) with 
the normal SN la population, the ratio of this value to 

30%, which is particularly im- 



the one in P 10 is 



portant at low-z, where the delayed component becomes 
dominant. This ratio sets a constraint on the fraction of 
low-s progenitors. 

5.4. Delay-time distributions 

Ideally, the rate should be able to help constrain the 
SN la delay-time distribution as follows: 



SNR Ia (t) 



v I SFR(i')DTD(i 
hp 



t')dt', (6) 



where t 1 is the age of the universe, tp the time at which 
the first stars formed (effectively zero for the SNLS), 
DTD (id) is the delay-time (td) distribution and v is the 
number of SNe la that explode per unit stellar mass 
formed. This number is related to the progenitor mass 
range and IMF, \1/(M), as well as the efficiency rj or frac- 
tion of white dwarfs that explode as SNe la, and the 
range of WD progenitor masses, and M^ D : 



& ^>{M)dM 



(7) 



Eq. |6] can set constraints on the progenitors through 
the delay-time distribution if the SFH is known. We 
assum e here the parameterization of [Hopkins fc Beacoml 
(12001 for a SalA IMF (h=0.7): 



c FH , (0.017 + 0.13z)fe 
SFH0C l + (z/3.3)5-3 ' 



(8) 



A simple parametrization power-law DTD oc i™, has 
been used by several authors to model the SN la 
rate, n ~ — 1 is obtained in the SN la r ate stu d- 
ies made by iTotani et all (|2008[ ) : IMaoz et all (j2010al lbh: 
iMaoz fe Badenesl (120101) and i s consistent with the DD 
scenario. iPritchet et"afl (|2008[) obtain a different power- 
law of n ~ —0.5 investigating the SN la rate as a func- 
tion of sSFR in the SNLS. It is interesting to see how the 
sub-sample of low-s behaves, in particular as the inferred 
DTDs tend to have a low tail at higher delay-times (after 
a sharp initial peak) where the subluminous SNe la are 
expected to lie. 

By fitting this model to the local and SNLS data up 
to z = 0.6, we obtain n = —0.93 ± 2.13 reflecting the 
large uncertainties in the rates. With the corrected z > 
0.6 rates, the power-law becomes n — —2.67 ± 6.98 (see 
Figure I15|) . This is accomplished by increasing the other 
free parameter, the scale v, by an order of magnitude 
such that many more SNe la happen at early times, in 
order to compensate for the steady rates at higher-z. 



We als o use the s i mple empirical models em- 
ployed in Strolgcr et al.l ((200 4) : an exponential func- 
tion motivated by the d ouble-degenerate scenario 
(|Tutukov fc Yungelsonl 119941 ) and a Gaussian func- 
tion motivated by a DT D with a specific delay-time 
(jDahlen fc Franssonl 119991 ) . The Gaussian distribution 
can be a proper approach for sub-groups of SNe la, such 
as the low-s one, but not for the entire population, which 
has been shown to be either bimodal or slowly declining 
in time. We use the f ormalism of the unim odal skew- 
normal DTD of eq.6 of iStrolger et~aT1 {2010) as a gener- 
alization of these models encompassing a wide variety of 
shapes and dist ributions. 

The model of IStrolger et al.l (|2010D has three indepen- 
dent variables, We obtain for the local and z < 0.6 rates: 
f = 6.4±3.8, u = 1.6±5.2 and a = -81.2 ±26400, which 
translate to a Gaussian behaviour of characteristic delay- 
time r = (5.11 ± 5.65) Gyr and width a = 0.94 ± 3.16. 
This fit is shown in Figure [15j It indicates that low- 
s SNe la come from progenitors 3-10 Gyr old following 
a single gaussian-shaped distribution. Such delay-times 
correspond to masses of 1.5 — 2.5M Q using t ~ m -2 5 . 
This sustains the argument that low-s SNe la come exclu- 
sively from a delayed channel of SNe la. Such estimates 
are higher t han t he lower delay-time limit obtained by 
iSchawinskH (|2009D based on individual SN la studies in 
early-type galaxies with opti cal and UV-data. O n the 
other hand, they agree with iBrandt et aD (|2010f ) who 
find that SDSS SNe la with s < 0.92 have delay-times of 
> 2.4Gyr. Fitting all points, corrected SNLS and low-z, 
we obtain r = (5.88 ± 4.84) Gyr and a = 0.15 ± 2.41 
(£ = 5.8 ± 3.6, lj = 0.2 ± 4.0 and a = 0.9 ± 17300) in- 
stead. This approximates to a much narrower Gaussian 
spanning 5.3-6.3 Gyr (1.95 — 2.10M©), time at which all 
low-s SNe la are created. This explains the sharp drop 
in the modeled rate at z ~ 0.9 seen in Figure [T5l Nev- 
ertheless, due to the large error uncertainties, all these 
results are poorly constrained. 

The explosion fraction is around 0.0001 — 0.0003 low- 
s SNe la per M Q for all models except for the power- 
law fit to all corrected data (t~ 2 ' 7 ) which is 0.001 low- 
s SNe la per M©. These fractions are lower than for 
the normal- s population -generally higher than ~ 10 ~ 3 
(e.g. iDahlen et al.l 12004 IMaoz et al.ll2010b[ ). Both pop- 
ulations differ by an order of magnitude which can be 
due either to a) a lower conversion efficiency rj of WD 
into SN la, or b) a lower progenitor mass range, tak- 
ing into account that low mass stars are more frequent, 
or a combination of both. If one assumes that the con- 
version efficiency is constant over mass -this can be the 
case fo r certain considerati ons of the SD and DD sce- 
narios (jPritchet et al.ll2008D - the allowed mass range of 
low-s SNe la would need to be extremely narrow. With 
a SalA IMF, and assuming that all main-sequence stars 
with 0.8M© < M < 8M Q become WDs, the cutoff value 
between low-s and normal-s progenitors (from Eq. [7]) 
would lie at ~ lAf Q , an extremely narrow range of 
stars with too long delay-times (> 10 Gyr) to be seen 
at z > 0.3. Changing the range of progen itors masses 
to 3M < M < 8M Q (jNomoto et al.l 119941 ) . the cutoff 
would need to be at ~ 3.2M©. On the other hand, us- 
ing the lower WD mass limit of 1.95M Q and cutoff at 
2.10M Q , from the Gaussian distribution previously ob- 
tained, we also get a consistent explosion fraction ratio 
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Fig. 15.— Best fit DTDs to the low-s SN la rate. Up to z < 0.6: 
power-law (dashed blue) with t~ 0-93±2.i3^ an( j Gaussian behaviour 
(solid purple) with r = (5.11 ± 5.65)Gyr and a = 0.94 ± 3.16. All 
corrected rates: power-law (dot-dashed blue) with j— 2 -67±6.98 ant j 
Gaussian behaviour (long dashed purple) with t = (5.88±4.84)Gyr 
and <r = 0.15 ± 2.41. The dotted gree n lines represent the best A 
model and errors for z < 0.6 (see fe.3j l, 

of low-s to normal- s. 

5.5. Extremely low-s SNe la 

We note that our sample does not consist of extremely 
subluminous SNc la. Looking at the local sample of typi- 
cal subluminous objects, such as SN1991bg or SN1999by, 
the stretch is on the low end, i. e., s = 0.5 — 0.65. In the 
SNLS 0.1 < z < 0.6 range, we only obtain one such 
object (see Figure [16]). This SN has a spectroscopic red- 
shift of z = 0.57 with LC parameters: s = 0.63 ± 0.08 
and c = 0.02 ± 0.06. The candidate is hosted in a galaxy 
of mass M = 10 WA7 M Q and sSFR = 10- 1193 Af o yr- 1 , 
typical of very low-stretch objects. 

Elfective day 

-20.0 0.0 20.0 




Fig. 16. — Extremely low-s candidate, 04D4ms, at spectroscopic 
redshift z = 0.57 with s = 0.63 ± 0.08 and c = 0.02 ± 0.06. 

This single SN corresponds to a single z = 0.35 rate of 
(0.25±°;^(stat)t°;™(MC)) x 10- 6 yr _1 Mpc- 3 ( see Fig- 
ure [17]). Despite the obvious uncertainty in this number, 



it is tempting to argue that perhaps this "extremely low- 
s'' population, i.e. s < 0.65, actually behaves differently 
than our s < 0.8 group. If we maintain the extremely 
low-s definition for th e local luminosity-function sam- 
ple of ILi et all (|2010al) as well, we obtain 5 objects with 
s < 0.65 (5.16 corrected for completeness), which corre- 
sponds to 6.9^^9% of their rate (if we were to use the 
entire local sample, 15 objects have s < 0.65, i.e., a ~ 9% 
fraction), that is a rate of (2.0 ± 1.1) x 10" 6 yr" 1 Mpc -3 . 
We note that this is much smaller than the 91bg frac- 
tions generally quoted, such as the 17.9^2% by ILi et ail 
(|2010al) , indicating that this photometric definition is 
more restrictive than their definition (which is also based 
on spectra). Contrarily, if we were to use their fraction 
of "subluminous" SNe la to estimate the stretch limit, 
we would obtain s ~ 0.75. We emphasize therefore the 
importance of having a clear and consistent definition 
of "subluminous" SNe la across the redshift range. As 
shown in Figure [T7l the rate evolution for this population 
also reveals a decline, perhaps even more rapid than for 
our standard low-s sample. 
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Fig. 17. — Single bin rate evolution for low-s candidates (s < 
0.8) and extremely low-s candidates with s < 0.65 from the SNLS 
(green and red circl es respectively), and the estimated fraction of 
the SN la rate from Li et al. (2010b) (triangles). Horizontal error 
bars indicate the width of the z-bin. 



6. SUMMARY 

We have conducted a search in the SNLS for low- 
stretch SNe la, defined here as objects with s < 0.8. 
Typical "subluminous" surveys considered more extreme 
objects based also on spectroscopic features which have 
equivalent stretches of < 0.72 — 0.75. The definition 
chosen suits better our constructed color-stretch and 
magnitude-color-stretch relations trained with low-z and 
SNLS objects. A lower transition value, nevertheless, 
provides also consistent LC relations. This argues for 
a continuous distribution of LC properties of SNe la, 
with low-s SNe la having distinctive redder colors and 
a steeper color-stretch relation. The exact stretch value 
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at which this transition occurs is not a sharp one and lies 
somewhere between and 0.69 — 0.87. 

The LC relations serve for the typing of low-s SNe la at 
high redshift and for modeling the survey efficiency. We 
find 18 photometrically-identified low-stretch SNe la up 
to z ~ 0.6, indicating that the primary reason for their 
previous non-detections at high redshifts was due to se- 
lection effects. Most of the candidates at z < 0.6 have 
s > 0.65 and colors similar to the local sample. Only one 
of them is comparable to the faintest local counterparts 
such as 1991bg. The candidates are observed to be pri- 
marily hosted in passive galaxies, with no star formation, 
as it is found locally. The spectroscopic candidates have 
no conclusive features due to low signal-to-noise. 

We also find a population of blue low-s SNe la at 
z > 0.6 hosted mostly in quiet host envirnoments. They 
represent the blue end of the low-s SN la sample at high- 
rcdshift, whereas its reddest counterpart is not seen due 
to selection effects. A 30% fraction of these have un- 
usual blue colors that are mostly attributed to the larger 
photometric uncertainties at high-z. Nonetheless, some 
of them occur in star-forming galaxies and are possible 
CC-SN contaminants (< 10%). Alternatively we specu- 
late that they could represent a different set of objects 
that evolve rapidly but have bluer colors than typical 
low-s SNe. 

We calculate a rate evolution of low-s SNe la as a 
function of redshift and the results indicate a steady or 
slightly decreasing trend up to z = 0.6. We correct the 
rates at z > 0.6, where our sample is incomplete, and find 
the evolution to be different from the increase of the nor- 
mal SN la population. The missclassification of CC-SNe 
as SNe la (or different unknown blue high-z transients) 
is an important source of systematics in the rate mea- 
surements. This would overestimate the rate and should 
become larger with redshift, so that the effect of a de- 
creasing rate would be even stronger. The main findings 
of this study are therefore not influenced by contamina- 
tion. 

An A-component model, proportional to just the mass 
of the galaxy, is sufficient to fit the rate. This is ex- 
pected as most of the candidates are hosted in passive 
galaxies with no star formation. Low-s SNe la are there- 
fore a population coming from evolved and/or metal- 
ri ch progenitors not d ependent on the SFH. As shown 
in iHowell et al.l ()2007l ) , this has consequences for cosmo- 
logical studies because the relative proportion of low to 
high stretch SNe la evolves at high redshift. The present 
study supports their findings at the extreme low-stretch 
end of the population. We find that A + B models have 
an A-component with ~ 30% coming from the low-s pop- 
ulation. 

Studies of the DTD of different scenarios help constrain 
progenitors masses and models, as well as the conversion 
efficiency of WDs int o SNe la. Despite th e large uncer- 
tainties, the model of iStrolger et al.l (|2010f) fitted to our 
data hints to a Gaussian distribution centered around 



r ~ 5 — 6 Gyr. The explosion fraction for those models is 
v ~ 10~ 3 — 10~ 4 low-s SNe la per Mq, an order of magni- 
tude lower than for normal-s objects. These results again 
argue for the low-s being part of a delayed component. 
Low-s SNe la, together with other long delay SNe la of 
higher-s, could come from a separate progenitor channel. 
Although both, the SD and DD, scenarios can produce 
delayed SNe la, the merger of two low-mass WDs produc- 
ing a sub-Chandrasekhar explosion seems more favorable 
to produce faint SNe in old stellar environments. 

Finally, future surveys should be able to determine the 
rate evolution of low-s much better and to higher red- 
shifts where the systematic errors in this study become 
more significant. Spectroscopic redshifts of high- z low- 
s SNe la are essential for this and could provide addi- 
tional hint on their properties and their evolution. Fur- 
ther comparison of different SN la populations will give 
clues to understand the formation and explosion mecha- 
nisms. 
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TABLE 1 

Nearby low-s objects that passed the observation and fit culls used in the analysis 



Name 


Redshift 


Stretch 


Source 


SN1986g 


0.002 


0.71 


Phillips ct al. (1987) 


SN1990af 


0.050 


0.75 


Hamuv et al. (1996) 


SN1991bg 


0.004 


0.50 


Turatto et al. (1996),FilipDenko et al. (1992) 


SN1992M 


0.044 


0.80 


Hamuv et al. (1996) 


SN1992bo 


0.019 


0.79 


Hamuv et al. (1996) 


SN1992br 


0.088 


0.60 


Hamuv et al. (1996) 


SN1993h 


0.024 


0.66 


Hamuv et al. (1996).Altavilla et al. (2004) 


SN1994m 


0.023 


0.76 


Riess et al. (1999).Altavilla et al. (2004) 


SN1995ak 


0.023 


0.79 


Riess et al. ( 1999) 


SN1998bp 


0.010 


0.63 


Jha et al. (2006) 


SN1998co 


0.018 


0.65 


Jha et al. (2006) 


SN1998de 


0.017 


0.64 


Modiaz et al. (2001) 


SN1999bm 


0.143 


0.68 


Kowalski et al. (2008) 


SN1999by 


0.002 


0.61 


Garnavich et al. (2004) 


SN1999da 


0.013 


0.59 


Krisciunas ct al. (2001) 


SN2000dk 


0.017 


0.77 


Jha et al. (2006) 


SN2001da 


0.017 


0.61 


Hicken et al. (2009) 


SN2002cx 


0.024 


0.57 


PhilliDS et al. (2007) 


SN2002hw 


0.018 


0.73 


Hicken et al. (2009) 


SN2003ic 


0.056 


0.77 


Hicken et al. (2009) 


SN2003iv 


0.034 


0.74 


Hicken et al. (2009) 


SN2005am 


0.008 


0.71 


Hicken et al. (2009). Li et al. (2006) 


SN2005M 


0.024 


0.60 


Taubenbereer et al. (2008) 


SN2005hf 


0.043 


0.76 


Hicken et al. (2009) 


SN2005ke 


0.005 


0.69 


Hicken et al. (2009) 


SN2005mc 


0.025 


0.69 


Hicken et al. (2009) 


SN2006ak 


0.038 


0.78 


Hicken et al. (2009) 


SN2006al 


0.068 


0.78 


Hicken et al. (2009) 


SN2006br 


0.025 


0.78 


Hicken et al. (2009) 


SN2006bw 


0.030 


0.71 


Hicken et al. (2009) 


SN2006bz 


0.028 


0.53 


Hicken et al. (2009) 


SN2006gj 


0.028 


0.71 


Hicken et al. (2009) 


SN2006gt 


0.045 


0.53 


Hicken et al. (2009) 


SN2006hb 


0.015 


0.69 


Hicken et al. (2009) 


SN2006je 


0.038 


0.48 


Hicken et al. (2009) 


SN2006mo 


0.037 


0.73 


Hicken et al. (2009) 


SN2006n 


0.014 


0.77 


Hicken et al. (2009) 


SN2006nz 


0.038 


0.62 


Hicken et al. (2009) 


SN2006ob 


0.059 


0.75 


Hicken et al. (2009) 


SN2006td 


0.016 


0.79 


Hicken et al. (2009) 


SN2007al 


0.012 


0.47 


Hicken et al. (2009) 


SN2007au 


0.021 


0.67 


Hicken et al. (2009) 


SN2007ax 


0.007 


0.53 


Hicken et al. (2009) 


SN2007ci 


0.018 


0.79 


Hicken ct al. (2009) 



TABLE 2 

Nearby non-Ia SNe with acceptable SN Ia SiFTO fits, enough data coverage and fitted s < 0.8. 



Name 


Type 


Redshift 


Stretch 


Source 


SN1999cu 


SNIIP 


0.004 


0.45 


Pastorello et al. (2004) 


SN2001b 


SNIb 


0.005 


0.60 


Tsvetkov (2006) 


SN2002ao 


SNIbn 


0.005 


0.54 


Pastorello et al. (2008) 


SN2002ap 


SNIc 


0.002 


0.68 


Gal- Yam et al. (2002); Folev et al. (2003); Yoshii et al. (1988) 


SN2003jd 


SNIc 


0.019 


0.73 


Valenti et al. (2008) 


SN2006aj 


SNIc 


0.033 


0.73 


Pian ct al. (2006); Sollcrman ct al. (2006) 
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TABLE 3 

Culls applied to the training sample 



Cull 




Number 






confirmed 


confirmed 




SNe Ia(low-s) 


CC-SNc(low-s) 




low- 2 


SNLS 


low- 2 


SNLS 


None 


246(57) 


112(3) 


37(12) 


47(9) 


observation 


207(44) 


106(3) 


14(6) 


21(4) 


XsiFTO an d sna kcs 


184(41) 


97(3) 


2(2) 


1(0) 



TABLE 4 
Culls applied to the SNLS sample 



Cull 






Number 






all 


low-s 


confirmed 


confirmed 








SNe la 


CC-SNc 


None 


5234 


1497 


407 


86 


observation 


2615 


837 


343 


62 


XsiFTO 


1054 


144 


338 


22 


- A -cstimatc_sn 


925 


118 


329 


8 


y? 

^ J,estimate_sn 


848 


109 


303 


5 


snakes 


710 


70 


295 


1 


FINAL 


710 


70 


295 


1 


2 < 0.6 


174 


18 


118 






TABLE 5 

Volumetric rates for the subluminous SNe Ia in the SNLS 
First and second column show the redshift at the middle of the bin and the mean redshift respectively. The third column 

IS THE NUMBER OF OBJECTS PER 2-BIN AND THE FOURTH IS THE MEAN WEIGHT (Wi = (1 + Zj)/(ti X ATj), WHERE €j is the efficiency) 
USED FOR THE POISSON ERRORS. FIFTH COLUMN INDICATE THE MEASURED RATE WITH STATISTICAL WEIGHTED POISSON ERRORS, AND SIXTH 
COLUMN THE MONTE-CARLO RATES WITH WEIGHTED POISSON, MONTE-CARLO AND CONTAMINATION ERRORS. THE SEVENTH COLUMN SHOWS 
THE CORRECTION FACTOR (WITH ERRORS) APPLIED TO GET THE LAST COLUMN HIGH-Z RATES CORRECTED FOR THE COLOR DISTRIBUTION AT 

2 < 0.6. 



^•mid 


< 2 > 


N 


< w >, a w 


^racas 

(10- 5 yr- 1 Mpc- 3 /i? n ) 


IMC 

(10- 5 yr- 1 Mpc- 3 ^ n ) 


Jcorr 


^corr 

(10- 5 yr- 1 Mpc- 3 ^ n ) 


0.20 


0.17 


2 


1.40, 0.75 


46 +U ' ti8 


0.42+; , ) ;^(stat)+; , ) ' , ^(MC) ± 0.08(cont) 






0.40 


0.43 


9 


1.08, 0.24 


o.74±°,;« 


°- 51 ±a2l( stat )±ao8( MC ) ± O-H(cont) 






0.60 
0.80 


0.63 
0.78 


19 
18 


1.16, 0.21 
1.44, 0.65 


0.6512:1? 
0.44t°;i; 


0-49l°; 2 °(stat)l°;J 2 (MC) ± 0.13(cont) 
0.40±g;^(stat)±°;^(MC) ± 0.12(cont) 


1.42 ±0.20 
1.89 ±0.26 


fiQ+ () - 35 
u - Dy -0 .29 

7^+ - 51 
u - '°-0.43 


0.95 


0.94 


11 


2.37, 3.02 


0-48t°; 3 2 ° 


0.27l[ ) ) ;g(stat)l° ) ;J®(MC) ± 0.08(cont) 


1.82 ±0.40 


49+0-62 



